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ABSTRACT 

River-based  hydrokinetic  turbine  power  generation  systems  have  been  studied  to 
introduce  an  effective  energy  flow  control  method.  Hydrokinetic  turbine  systems  share  a 
lot  of  similarities  with  wind  turbine  systems  in  terms  of  physical  principles  of  operation, 
electrical  hardware,  and  variable  speed  capability  for  optimal  energy  extraction.  A  multi¬ 
pole  permanent  magnet  synchronous  generator  is  used  to  generate  electric  power  because 
of  its  ability  to  reach  high  power  density  and  high  thrust  at  low  speed.  A  3-phase  diode 
rectifier  is  used  to  convert  AC  power  from  the  generator  into  DC  power  and  a  boost 
converter  is  used  to  implement  energy  flow  control.  On  the  load  side,  an  electronic 
voltage  load  is  used  for  test  purposes  to  simulate  a  constant  DC  bus  voltage  load,  such  as 
a  battery.  A  dynamic  model  of  the  entire  system  is  developed  and  used  to  analyze  the 
interaction  between  the  mechanical  structure  of  water  turbine  and  electrical  load  of  the 
system,  based  on  which  a  maximum  power  point  tracking  control  algorithm  is  developed 
and  implemented  in  the  boost  converter.  Simulation  and  experimental  results  are 
presented  to  validate  the  proposed  MPPT  control  strategy  for  hydrokinetic  turbine 
system.  Similar  to  the  wind  turbine  system,  hydrokinetic  turbine  system  usually  requires 
a  gear  box  to  couple  the  turbine  and  the  generator  because  the  operating  speed  range  for 
the  hydrokinetic  turbine  is  much  lower  than  the  operating  speed  range  for  most  PMSGs. 
However,  the  gear  box  coupling  adds  additional  transmission  power  losses.  Therefore  a 
high-thrust  low-speed  permanent  magnet  synchronous  generator  is  designed  to  couple 
with  the  water  turbine  without  a  gear  box. 
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1.  INTRODUCTION 


1.1.  MOTIVATION 

As  a  result  of  rising  concerns  over  greenhouse  gas  emissions,  the  environment, 
rising  electricity  demand  and  exhaustion  of  fossil  energy  sources,  trends  in  energy 
production  have  become  more  and  more  focused  on  using  renewable  energy  resources 
such  as  wind,  solar,  hydro,  biomass  and  geothermal  [1-3].  Among  these  renewable 
energy  resources,  hydro  and  wind  power  have  great  potentials  of  development  to  be  able 
to  compete  with  the  classical  fossil  energy  sources  such  as  coal,  oil  and  natural  gas  [4] . 
Wind  power  over  the  past  decade,  has  undergone  a  substantial  advancement  in  terms  of 
control  strategies,  custom  designed  generators,  maintenance  cost  and  cost 
competitiveness,  becoming  more  and  more  important  in  renewable  energy  market  [5]. 
Hydroelectricity  has  been  utilized  all  over  the  world  for  decades,  usually  in  the  form  of 
high  power  hydroelectric  stations  which  require  the  construction  of  impoundments  or 
diversion  of  conventional  hydroelectric  facilities  based  on  dams  or  penstocks.  Nowadays, 
since  the  suitable  sites  for  these  plants  are  rare  and  concerns  over  the  effects  of  these 
plants  on  local  environment  rise,  it  is  no  more  possible  to  set  up  such  plants  in  most 
European  and  North-American  countries.  Nonetheless  substantial  capacities  of  these 
hydroelectricity  plants  remain  in  Africa,  Asia,  and  South  America  [6].  In  order  to  utilize 
the  vast  potential  of  hydro  energies  without  the  construction  of  dams  and  penstocks, 
hydrokinetic  turbine  energy  conversion  system  (HTECS)  becomes  more  and  more 
popular  and  is  at  the  forefront  of  pre-commercial  or  commercial  deployment  at  various 
locations  in  the  world  [7]. 

Hydrokinetic  turbines  share  many  similarities  with  wind  turbines  in  terms  of 
physical  principles  of  operation,  electrical  hardware,  and  variable  speed  capability  for 
optimal  energy  extraction  [8].  Therefore,  HTECS  can  benefit  in  many  ways  from  the 
many  years  of  advanced  development  in  the  wind  turbine  energy  conversion  system 
(WTECS).  However,  significant  differences  between  HTECS  and  WTECS  do  exist,  such 
as  a  narrower  tip  speed  ratio  range,  cavitation  limits  and  the  requirement  to  deal  with 
harsh  environments  for  HTECS  [9].  One  of  the  most  significant  differences  between 
WTECS  and  HETCS  is  that  wind  speed  varies  dramatically  over  timescales  of  seconds  to 


2 


minutes  while  water  speed  tends  to  be  relatively  constant  with  seasonal  variations.  Since 
the  maximum  output  power  point  for  the  same  wind  turbine  varies  under  different  wind 
speeds,  the  maximum  power  point  tracking  (MPPT)  control  for  wind  turbines  has  been 
the  subject  of  detailed  research  [10-12].  Though  water  speed  tends  to  be  relatively 
constant  with  small  seasonal  variations,  the  fact  that  water  density  is  more  than  800  times 
bigger  than  air  density  means  that  even  a  small  water  speed  variation  would  yield  a 
substantial  output  power  change.  Therefore,  a  maximum  power  point  tracking  strategy 
for  the  hydrokinetic  turbine  system  is  necessary  to  be  able  to  extract  the  maximum 
amount  of  potential  power  from  the  hydrokinetic  turbine. 

In  order  to  convert  the  mechanical  hydrokinetic  energy  into  electrical  power, 
electrical  generators,  such  as  permanent  magnet,  synchronous  and  induction  generators, 
are  used.  For  small  scale  HTECS  and  WTECS,  permanent  magnet  synchronous 
generators  (PMSG)  are  widely  used  because  of  their  high  reliability,  simple  structure,  low 
noise  and  high  power  density  [13].  Then  the  electrical  power  is  processed  using  power 
electronic  devices  such  as  ac-dc  rectifier  (controlled/uncontrolled),  dc-dc  converter  and 
dc-ac  inverter  [14].  Based  on  whether  the  system  is  interfaced  with  grid  or  not, 
hydrokinetic  turbine  systems,  like  wind  turbine  systems,  have  two  types  of  configurations 
in  general.  One  is  grid-tie  hydrokinetic  turbine  system  and  the  other  one  is  stand-alone 
hydrokinetic  turbine  system  as  shown  in  Fig.  1.1  and  Fig.  1.2  respectively.  For  grid-tie 
hydrokinetic  turbine  systems,  a  dc-ac  inverter  is  necessary  to  convert  dc  power  into  ac 
power  to  interface  with  the  grid.  For  stand-alone  hydrokinetic  turbine  systems,  battery 
banks  are  usually  used  to  store  the  extracted  energy  as  well  as  to  provide  a  constant  dc 
bus  voltage  on  the  load  side.  DC  battery  bank  can  store  extra  energy  when  the  load  power 
demand  is  less  than  the  available  power  from  the  hydrokinetic  turbine  and  output  power 
when  the  load  power  demand  exceeds  the  available  power  from  the  hydrokinetic  turbine 
[13,  15].  Stand-alone  hydrokinetic  systems  are  common  in  rural  areas  where  there  are  no 
accesses  to  a  power  grid  [13]. 
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Figure  1.1.  Grid-tie  hydrokinetic  turbine  system 
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Figure  1.2.  Stand-alone  hydrokinetic  turbine  system 


1.2.  PROPOSED  CONTROL  APPROACH 

Generally,  the  objective  of  controlling  a  turbine  system  is  to  maximize  the  energy 
yield  because  with  different  water  (or  wind)  velocity,  the  maximum  power  point  for  the 
hydrokinetic  (or  wind)  turbine  occurs  at  different  turbine  rotor  speed  [16].  In  medium  and 
large  turbine  systems,  pitch  angle  control  is  usually  adopted  to  optimize  output  power  as 
well  as  to  reduce  torque  and  output  power  variations  for  high  water  (or  wind)  speed  [17]. 
However,  for  small  size  turbine  systems,  the  mechanical  structure  limits  make  pitch  angle 
control  unrealistic  [13].  Therefore,  a  power  converter  is  usually  adopted  to  implement  the 
MPPT  control  in  small-scale  turbine  systems,  such  as  boost  converter.  In  order  to  control 
the  operation  of  the  power  converter,  a  microcontroller  is  needed. 
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Conventionally,  in  order  for  the  turbine  to  operate  at  maximum  power  point,  the 
turbine  power  versus  rotational  speed  characteristic  for  a  specific  turbine  is  stored  in  a 
memory.  By  measuring  the  rotational  speed  of  the  turbine,  the  difference  between  the 
maximum  power  available  and  the  actually  output  power  can  be  found.  This  error 
between  the  two  power  values  then  feeds  back  to  the  control  of  the  generator  so  that  it 
can  operate  at  the  speed  that  will  generate  maximum  power  [13].  However,  pre¬ 
measuring  the  turbine  power  versus  speed  characteristic  is  costly  for  small-scale  turbine 
systems.  Therefore,  a  dynamic  MPPT  control  strategy  for  a  small-scale  hydrokinetic 
turbine  system  is  proposed  in  this  thesis.  This  proposed  MPPT  control  approach  enables 
the  system  to  dynamically  track  the  maximum  power  point  of  the  hydrokinetic  turbine 
without  pre-measuring  the  power  versus  rotational  speed  characteristics  of  the 
hydrokinetic  turbine.  Figure  1 .3  shows  the  system  configuration  of  the  stand-alone 
hydrokinetic  turbine  system  with  the  proposed  MPPT  control  implemented  in  a  boost 
converter. 
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Figure  1.3.  Stand-alone  hydrokinetic  turbine  system  with  MPPT  control  implemented  in 

the  boost  converter 


5 


In  this  control  method,  the  output  voltage  and  current  are  constantly  measured  and 
updated  by  the  microcontroller,  which  would  give  the  value  of  the  output  power.  Based 
on  the  change  of  the  output  power,  the  duty  ratio  will  be  changed  accordingly.  Since  the 
output  voltage  is  constant,  varying  the  duty  ratio  of  the  boost  converter  will  cause  the 
output  voltage  of  the  rectifier  to  change  which,  as  a  result,  will  change  the  rotational 
speed  of  the  generator. 

1.3.  DOCUMENT  ORGANIZATION 

This  thesis  is  organized  into  four  sections.  Section  1  provides  a  general 
introduction  for  the  hydrokinetic  system  and  a  proposed  MPPT  control  approach  for 
hydrokinetic  turbines. 

Section  2  describes  the  system  level  design  of  hydrokinetic  energy  conversion 
systems.  It  covers  different  methods  and  system  configurations  for  hydrokinetic  energy 
conversion  systems  such  as  vertical  and  horizontal  turbine  systems,  ducted  and  unducted 
turbine  systems  as  well  as  vortex  induced  vibration  systems. 

Section  3  presents  a  dynamic  model  of  the  hydrokinetic  turbine  system  based  on 
which  a  MPPT  control  method  is  developed.  Simulation  in  Matlab1  and  experimental 
results  are  also  presented  in  this  section  to  validate  the  proposed  control  method 

Section  4  describes  a  high-thrust  low-speed  permanent  magnet  synchronous 
generator  design  for  the  hydrokinetic  turbine  systems 

Section  5  presents  the  conclusions.  Future  work  and  project  extensions  for  the 
hydrokinetic  system  are  also  discussed  in  this  section 

Appendix  A  contains  the  schematic  and  board  layout  for  the  3  phase  rectifier  and 
boost  converter  used  for  the  experiment  data  collection. 

Appendix  B  contains  the  C  codes  of  the  developed  MPPT  algorithm  implemented 
in  the  microcontroller. 


1  Matlab  is  a  registered  trademark  of  The  MathWorks,  Inc. 
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2.  SYSTEM  LEVEL  DESIGN  OF  HYDROKINETIC  TURBINE  SYSTEMS 

2.1.  HYDROKINETIC  ENERGY  UTILIZATION 

Hydrokinetic  energy  generation  technologies  have  been  under  development  for 
decades,  but  this  development  has  been  sporadic  and  inconsistent.  Traditionally, 
hydrokinetic  energy  can  be  categorized  as  wave  energy  and  current  energy  such  as  tidal 
current,  open-ocean  current  and  river  current.  Wave  energy  technologies  are  mostly  used 
in  marine  hydrokinetic  systems.  Wave  energy  converters  that  are  currently  under 
deployment  can  be  categorized  as  absorbers,  attenuators,  oscillating  water  columns, 
overtopping  devices,  and  inverted  pendulums  [18].  Hydrokinetic  current  energy,  which  is 
common  in  both  oceans  and  rivers,  is  usually  converted  into  electrical  energy  by 
hydrokinetic  turbines  just  like  the  wind  kinetic  energy  is  converted  into  electrical  energy 
by  wind  turbines.  In  addition,  hydrokinetic  current  energy  can  be  converted  into  electrical 
energy  by  Vortex  Induced  Vibration  (VIV)  systems,  in  which  a  cylinder  that  is 
submerged  into  the  water  would  vibrate  as  the  water  flow  hits  the  cylinder  and  this 
vibration  motion  is  converted  into  electricity  by  a  linear  generator. 

2.2.  HYDROKINETIC  TURBINE  SYSTEM  DESIGN 

2.2.1.  Horizontal  and  Vertical  Hydrokinetic  Turbine  System.  The  first  step  to 
design  a  hydrokinetic  turbine  system  (HTS)  is  to  decide  whether  to  use  horizontal  turbine 
system  or  vertical  turbine  system.  For  vertical  hydrokinetic  turbine  system,  the  biggest 
advantage  is  that  it  can  have  its  gearbox  and  generator  above  water  level.  However,  based 
on  the  experiments  conducted  in  [19],  the  fact  that  vertical  hydrokinetic  turbine  system 
can  have  its  gearbox  and  generator  above  water  level  does  not  seem  to  have  any 
advantages  over  the  horizontal  hydrokinetic  turbine  because  the  gearbox  and  generator  in 
the  vertical  turbine  system  would  have  to  be  as  watertight  as  the  horizontal  turbine 
system.  For  horizontal  hydrokinetic  turbine  system,  the  biggest  advantage  is  that  it  is  a 
conventional  turbine  system  and  its  characteristics  have  been  studied  for  decades 
especially  in  the  wind  power  industry.  In  addition,  it  is  easier  for  horizontal  turbine 
systems  to  integrate  a  water  speed  increaser  to  increase  the  water  velocity  interacting 
with  the  turbine  blades  whereas  for  vertical  hydrokinetic  turbine  systems,  the  water 
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velocity  increase  may  have  to  be  directed  at  a  proper  angle  based  on  the  turbine  blade 
design  to  improve  the  power  output  [19].  Figure  2.1  shows  a  ducted  horizontal 
hydrokinetic  turbine.  With  the  help  of  duction  to  increase  the  water  velocity,  the  need  for 
a  gearbox  might  be  eliminated. 


Figure  2.1.  Lunar  energy’s  rotech  tidal  turbine  (taken  from  LunarEnergy-URL  2007) 


2.2.2.  Permanent  Magnet  Generator  and  Induction  Generator  System.  While 
hydrokinetic  turbine  converts  the  kinetic  energy  stored  in  the  water  mass  into  mechanical 
energy  available  at  the  end  of  turbine  shaft,  a  device  is  needed  to  convert  this  mechanical 
energy  into  electrical  energy  for  further  use.  Generally,  two  kinds  of  generators  are  used 
in  wind  and  hydrokinetic  turbine  systems,  permanent  magnet  synchronous  generator  and 
induction  generators  such  as  doubly-fed  induction  generators  and  squirrel  cage  induction 
generators  [20].  Permanent  magnet  synchronous  generators  are  usually  employed  for  low 
power  systems  while  induction  machines  are  usually  employed  for  high  power  systems. 
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A  doubly-fed  induction  generator  (DFIG)  is  also  able  to  supply  power  at  constant  voltage 
and  constant  frequency  while  the  rotor  speed  varies,  which  makes  it  suitable  for  grid-tie 
variable-speed  turbine  systems  [21].  In  addition,  another  big  advantage  of  the  DFIG- 
based  system  is  that  the  power  electronic  equipment  only  needs  to  handle  a  fraction  of  the 
total  system  power,  usually  20-30%,  which  reduces  the  power  loss  and  the  cost  of  power 
electronic  devices  [22] . 

River-based  hydrokinetic  turbine  systems  are  usually  KW  level  systems  which  are 
categorized  into  low  power  turbine  systems.  Therefore  PMSG  is  usually  adopted  for 
electrical  power  generation  for  river-based  hydrokinetic  turbine  systems. 

2.2.3.  Power  Converters  and  Control  for  Hydrokinetic  Turbine  System. 
Power  converters  are  necessary  in  a  hydrokinetic  turbine  system  to  process  the  electrical 
power  generated  by  the  generator  and  in  the  mean  time,  control  the  hydrokinetic  turbine 
operations.  Depending  on  the  load  type  of  the  system,  different  power  converters  are 
needed. 

For  a  grid-tie  hydrokinetic  turbine  system  shown  in  Fig  1.1,  an  ac/dc  rectifier  is 
used  to  convert  the  ac  power  coming  out  of  the  generator  in  to  dc  power.  Then  a  dc/dc 
converter  is  connected  after  the  rectifier.  To  interface  with  the  grid,  a  dc/ac  inverter  is 
connected  in  between  the  dc/dc  converter  and  the  grid.  There  are  two  control  objectives 
for  grid-tie  hydrokinetic  turbine  system.  One  is  to  ensure  maximum  output  power  capture 
of  the  hydrokinetic  turbine.  The  other  one  is  to  ensure  the  power  flow  injecting  into  the 
grid  which  requires  a  constant  dc  bus  voltage  on  the  output  side  of  the  dc/dc  converter 
and  the  grid  side  current  to  synchronize  with  the  grid  voltage.  Since  there  are  two  power 
converters  in  this  kind  of  system  configuration,  therefore  there  are  two  different  control 
strategies.  The  first  control  strategy  is  illustrated  in  Fig  2.2.  The  turbine  output  power 
tracking  is  realized  by  the  dc/dc  boost  converter  while  the  dc  bus  regulation  and  grid 
current  synchronization  control  is  realized  by  the  dc/ac  inverter.  The  second  control 
strategy  is  illustrated  in  Fig  2.3.  In  the  second  control  strategy,  the  boost  converter  and 
the  dc/ac  inverter  have  shifted  functions  in  comparison  with  the  first  control  strategy.  The 
boost  converter  is  used  to  regulate  the  dc  bus  voltage  and  the  dc/ac  inverter  is  used  to 
track  the  turbine  output  power  and  synchronize  the  grid  current.  These  two  control 
strategies  for  PMSG  grid-tie  turbine  system  are  elaborated  in  [2],  and  the  simulation 


results  show  that  both  control  strategies  have  the  same  power  control  capabilities  under 
rapid  wind  or  water  speed  changes.  However,  the  second  control  strategy  has  a  better 
power  control  capability  compared  with  the  first  control  strategy  under  grid  fault. 
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Figure  2.2.  Control  strategy  I  of  grid-tie  hydrokinetic  turbine  system 
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Figure  2.3.  Control  strategy  n  of  grid-tie  hydrokinetic  turbine  system 


For  a  stand-alone  hydrokinetic  turbine  system  shown  in  Fig  1.2,  an  ac/dc  rectifier 
is  used  to  convert  the  ac  power  coming  out  of  the  generator  in  to  dc  power.  Then  a  dc/dc 
converter  is  connected  in  between  a  battery  load  and  the  ac/dc  rectifier.  This 
configuration  is  similar  to  the  grid-tie  system  except  for  that  a  dc/ac  inverter  is  no  longer 
needed  and  an  energy  storage  device,  such  as  battery  is  used.  Since  there  is  no  grid 
connection,  the  only  control  objective  for  a  stand-alone  hydrokinetic  turbine  system  is  to 
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ensure  the  maximum  power  capture  of  the  hydrokinetic  turbine  which  is  shown  in  Fig 
1.3. 

2.2.4.  System  Design  for  Missouri  S&T  Hydrokinetic  Project.  The  objective 
for  the  Missouri  S&T  hydrokinetic  project  is  to  design  a  prototype  of  a  river-based 
hydrokinetic  turbine  system.  The  estimated  output  power  of  the  system  is  7.7  w  and  the 
system  is  designed  to  produce  power  for  individual  household  without  grid  interface.  The 
output  power  is  limited  by  the  size  of  the  water  tunnel  and  the  hydrokinetic  turbine  that 
will  be  tested  in  the  water  tunnel.  Based  on  these  preliminary  designing  objectives,  a 
system  level  design  for  this  hydrokinetic  turbine  system  is  carried  out. 

First  of  all,  the  general  turbine  type  of  the  system,  vertical  or  horizontal,  has  to 
be  determined.  Even  though,  vertical  turbine  system  is  innovative  and  allows  the 
generator  to  be  mounted  above  the  water  level,  the  lack  of  predictability  of  its 
performance  and  the  difficulty  to  integrate  water  velocity  increaser  into  the  system  leads 
to  a  horizontal  axis  turbine  system.  Also,  because  of  the  difficulty  to  connect  the 
generator  to  the  turbine  in  the  water,  the  generator  is  mounted  on  top  of  a  platform  that  is 
placed  above  the  water  level.  This  requires  a  beveling  gear  to  convert  the  horizontal 
rotational  movement  into  vertical  rotational  movement. 

Once  the  turbine  type  is  determined,  a  generator  needs  to  be  selected  to  convert 
the  mechanical  energy  into  electrical  energy.  As  discussed  in  section  2.2.2,  permanent 
magnet  synchronous  generators  are  usually  employed  for  low  power  systems,  because  of 
its  high  reliability,  simple  structure,  low  noise  and  high  power  density.  In  addition,  for 
PMSG  to  operate  there  is  no  need  to  add  extra  field  excitation  winding  like  induction 
generators.  Therefore  PMSG  is  used  for  this  system.  The  selected  PMSG  needs  to 
coordinate  well  with  the  hydrokinetic  turbine  to  achieve  good  performance  of  the  overall 
system.  Generally  speaking,  there  are  two  important  parameters  of  the  turbine  and  the 
generator  that  needs  to  be  properly  coordinated.  One  parameter  is  the  turbine  rotational 
speed  and  the  generator  rotational  speed.  If  a  gearbox  is  not  used  in  between  the  turbine 
and  the  generator,  the  turbine  maximum  rotational  speed  should  not  exceed  the  generator 
maximum  rotational  speed.  In  addition,  the  rated  generator  speed  needs  to  be  close  to  the 
operating  speed  range  of  the  turbine.  For  example,  if  the  turbine  operating  speed  range  is 
10  times  less  than  the  rated  speed  of  the  generator,  the  overall  system  is  highly  inefficient 
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because  the  operating  speed  is  far  less  than  the  rated  speed  of  the  generator.  One  possible 
solution  is  to  use  a  gearbox  (speed  increaser)  in  between  the  turbine  and  the  generator,  so 
that  the  rotational  speed  for  the  generator  can  be  increased  to  its  rated  operating  speed. 
The  other  parameter  is  the  maximum  torque  of  the  turbine  and  the  starting  torque  of  the 
generator.  The  maximum  torque  of  the  turbine  has  to  be  greater  than  the  starting  torque  of 
the  generator  otherwise  the  turbine  is  not  able  to  drive  the  generator.  A  SHINANO 
KENSHI  PMSG  which  is  rated  at  80  W  and  3000  rpm  is  selected  for  this  system.  Since 
the  rated  speed  of  the  selected  PMSG  far  exceeds  the  operating  speed  range  of  the 
hydrokinetic  turbine,  a  1: 10  ratio  speed  increaser  is  built  for  the  system. 

After  the  generator  is  selected,  the  next  step  is  to  design  the  power  converters  for 
the  system.  As  shown  in  Fig.  1.2,  for  stand-alone  hydrokinetic  turbine  systems,  a  diode 
rectifier  and  a  boost  converter  is  usually  used  to  process  the  electrical  power  flow  and  a 
battery  is  needed  to  store  the  generated  electrical  energy. 

In  terms  of  control  for  the  stand-alone  prototype  hydrokinetic  turbine  system, 
maximum  power  point  tracking  control  is  adopted  to  ensure  the  maximum  capture  of  the 
potential  hydrokinetic  energy  which  makes  the  system  more  cost-effective.  The 
developed  MPPT  control  strategy  for  the  hydrokinetic  turbine  system  is  elaborated  in 
section  3. 
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3.  MPPT  CONTROL  OF  STAND-ALONE  HYDROKINETIC  TURBINE 

3.1.  INTRODUCTION 

A  hydrokinetic  turbine,  like  a  wind  turbine,  has  the  variable  speed  capability  for 
optimal  energy  extraction.  For  different  water  velocity,  the  maximum  power  point  for  the 
hydrokinetic  turbine  occurs  at  different  turbine  rotor  speed  [16].  Therefore,  in  order  to 
capture  the  maximum  amount  of  potential  energy  from  the  hydrokinetic  turbine,  the 
turbine  rotational  speed  has  to  be  optimized. 

In  medium  and  large  turbine  systems,  pitch  angle  control  is  usually  adopted  to 
optimize  output  power  as  well  as  to  reduce  torque  and  output  power  variations  for  high 
water  speed  [17].  However,  for  small  size  turbine  systems,  the  mechanical  structure 
limits  make  pitch  angle  control  unrealistic  [13].  Therefore,  a  power  converter  is  usually 
adopted  to  implement  the  MPPT  control  in  small-scale  turbine  systems,  such  as  boost 
converter. 

In  order  to  analyze  the  mechanical  and  electrical  interaction  of  the  hydrokinetic 
turbine  system,  detailed  mathematical  model  of  each  component  of  the  system  is  studied, 
based  on  which  an  MPPT  control  method  is  developed  in  this  section. 

3.2.  HYDROKINETIC  TURBINE  MODEL 

The  developed  model  is  based  on  the  steady-state  power  characteristics  of  the 
hydrokinetic  turbine.  The  output  power  of  a  hydrokinetic  turbine  is  given  by 

Pm=0.5pxrVCp(A,/])  (1) 

where  Pm  is  the  turbine  output  power,  p  is  the  water  density(997  kg/m3  at  25  °C), 
r  is  the  blade  radius(m),  v  is  the  water  velocity(m/s),  Cp  is  the  power  coefficient  of  the 
hydrokinetic  turbine,  X  is  the  tip  speed  ratio  and  p  is  the  pitch  angle  of  the  turbine. 

Therefore,  if  the  water  density  p ,  blade  radius  r  and  water  velocity  v  are  constant, 
the  output  power  Pm  is  proportional  to  the  power  coefficient  Cp  which  is  determined  by 
the  value  of  X  and  p.  However,  for  a  fixed  pitch  angle  hydrokinetic  turbine  which  is  the 


main  interest  of  discussion  in  this  thesis,  the  power  coefficient  is  only  determined  by  2 
Therefore,  the  output  power  of  a  hydrokinetic  turbine  is  given  by 
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Pm=0.5pxrVCp(V  (2) 

The  tip  speed  ratio  of  a  hydrokinetic  turbine  is  defined  as 


v 

where  r  is  the  blade  radius(m),  v  is  the  water  velocity(m/s),  and  com  is  the 
rotational  speed  of  the  hydrokinetic  turbine(rad/s). 

Based  on  the  hydrodynamic  analysis  of  the  hydrokinetic  turbine  used  in  this 
hydrokinetic  project  [26],  a  simplified  relationship  between  Cp  and  2  when  the  pitch  angle 
is  12°  is  given  by 
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where  B  is  the  number  of  blades  of  the  hydrokinetic  turbine,  which  is  3,  and  C  is 
the  ratio  of  lift  coefficient,  C/,  to  drag  coefficient,  Cci,  which  is  30.  All  of  the  turbine 
parameters  are  based  on  the  hydrodynamic  analysis  of  the  hydrokinetic  turbine  that  will 
be  used  in  this  project. 

Based  on  Eqn.  (4),  the  Cp-X  characteristic  of  the  hydrokinetic  turbine  is  shown  in 
Fig.  3.1.  The  maximum  value  of  Cp  is  0.395  which  is  achieved  at  2=5.2. 

Figure  3.2  gives  the  characteristics  of  the  hydrokinetic  turbine  power  coefficient 
Cp  for  different  values  of  water  speed  where  Vl=l  m/s,  V2=0.9  m/s,  V3=0.8  m/s,  and 
V4=0.7  m/s.  The  maximum  turbine  output  power  Pm_maximum  under  different  water 
velocities  is  achieved  at  different  rotational  speed  as  shown  in  Fig  3.3. 
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Figure  3.1.  Cp-X  characteristic  of  the  hydrokinetic  turbine 
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Figure  3.2.  Cp-com  characteristic  of  the  hydrokinetic  turbine  at  different  water  speed 
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Figure  3.3.  Pm-com  characteristic  of  the  hydrokinetic  turbine  at  different  water  speed 

3.3.  PMSG  AND  POWER  CONVERTER  CIRCUIT  MODEL 

As  shown  in  Fig.  1.2,  a  typical  stand-alone  hydrokinetic  turbine  system  consists 
of  a  permanent  magnet  synchronous  generator  which  converts  the  mechanical  power  into 
electrical  power,  a  diode  rectifier  which  converts  the  ac  power  into  dc  power,  a  dc-dc 
boost  converter  which  implements  the  control  and  battery  banks  which  store  the  energy 
generated.  In  order  to  analyze  the  interaction  between  the  mechanical  hydrokinetic 
turbine  and  the  electrical  part  of  the  system,  a  modeling  of  the  PMSG,  diode  rectifier  and 
the  boost  converter  is  necessary.  The  connection  schematic  of  the  PMSG,  diode  rectifier, 
boost  converter  and  the  battery  load  is  shown  in  Fig.  3.4.  Each  phase  of  the  PMSG  is 
modeled  as  back  electromotive  force  (EMF)  E,  stator  resistance  R  and  stator  inductance 
L.  Since  the  stator  resistance  R  is  usually  small,  it  is  neglected  for  the  analysis  in  this 


thesis. 
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Figure  3.4.  Connection  circuit  of  PMSG,  diode  rectifier,  boost  converter  and  battery  load 

The  generated  back  electromotive  force  E  of  PMSG  is  expressed  as 

E  =  Ka>m  (5) 

where  K  is  the  back  EMF  constant  of  the  PMSG,  com  is  the  rotational  speed  of  the 

PMSG. 

For  each  phase  of  the  PMSG,  the  phase  voltage  can  be  written  as 

E  =  jcoeLI  +  V  (6) 

In  terms  of  magnitude,  Eqn.  (6)  becomes 

E2  =  (coeLI)2  +V2  (7) 

The  relationship  between  the  mechanical  frequency  co,„  and  the  electrical 
frequency  coe  is  given  by 


=  npo)m 


(8) 
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where  np  is  the  number  of  pairs  of  poles  of  the  PMSG. 

Substituting  (5)  and  (8)  into  (7),  I  can  be  expressed  as 

I_-J(K‘»„f-V2  (9 

Lnpa,n 

Assuming  the  power  from  the  generator  is  converted  into  dc  power  through  diode 
rectifier  with  a  unity  power  factor  and  that  the  load  current  is  continuous  with  no  power 
losses,  this  yield 


Pg=3VI  =  VdId 


(10) 


The  mean  value  of  V,j  with  respect  to  Vu^,eak  is  given  by 


Vd  =-\\VLL_PeakC°Sdd0 


3 

V  =—  V 

d  LL  peak 

71 


(11) 

(12) 


where  VLL_peak  is  the  line  to  line  peak  voltage  of  the  PMSG. 

Therefore  the  relationship  between  V,j  and  Vll  and  V  can  be  expressed  as 


vd=-jivu 


n 


K  =-Sv 


n 


(13) 

(14) 


where  Vll  is  the  rms  value  of  the  line  to  line  voltage  of  the  PMSG  and  V  is  the 
line  to  neutral  voltage  of  the  PMSG. 
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Substituting  (14)  into  (10),  the  relationship  between  ld  and  1  can  be  expressed  as 


1  =  —  I 

d  s 


(15) 


From  (9),  (10),  (14)  and  (15),  Pgc an  be  expressed  as 


P„  = 


nV, 


SLnpa)m 


(Wj 


f  nVA  v 


v3>/6y 


(16) 


From  (16),  the  condition  for  which  the  PMSG  generates  power  is  obtained 


f  7  S 


>0 


(17) 


Solving  (17),  the  minimum  speed  for  which  the  PMSG  generates  power  is 


.  nVd 
m  > 


iSk 


(18) 


/t,  v  t 

Therefore,  the  rotational  speed  of  the  generator  has  to  be  higher  than  — j=—  for 
the  generator  to  generate  power. 

Figure  3.5  shows  the  PMSG  output  power  with  respect  to  the  generator  rotational 
speed  at  various  values  of  Vd.  The  minimum  generator  rotational  speed  for  which  the 
PMSG  starts  to  generate  power  is  determined  by  the  value  of  Vd .  The  higher  the  Vd ,  the 
higher  the  generator  speed  is  before  it  starts  to  export  energy.  Therefore,  in  order  to 
produce  energy  at  low  water  speed,  it  is  necessary  to  impose  a  lower  value  of  Vd  which 
can  be  achieved  by  using  a  boost  converter. 

Figure  3.6  shows  the  PMSG  output  power  with  respect  to  Vd  at  different  generator 
rotational  speed.  For  each  rotational  speed,  there  is  only  one  optimal  value  of  Vd  that 
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generates  the  maximum  power.  Therefore,  in  order  to  track  the  maximum  power  point, 
optimal  Vd  has  to  be  found. 

Since  a  boost  converter  is  used  in  between  the  diode  rectifier  and  the  load  side, 
the  following  Eqn.  can  be  obtained 


1 

K  1  ~D 


(19) 


where  D  is  the  duty  ratio  of  the  boost  converter. 
Substituting  Error!  Reference  source  not  found,  into 
Error!  Reference  source  not  found.,  Pg  can  be  expressed  as 


P  = 


x(l-D)Vk 


bus 


SLnpcom 


(K*J- 


n{\  -  D)V, 


bus 


3^6 


(20) 


450 


400 


.150  - 


;  300 

5 

S 

e 

o-  250 


200 


150 


100 


50 


/  j 

\A 


i 

IXj 


/[] 


/ 


u.  I 
20 


60  80  IOO 

Gcncralor  Rotational  Speed  <d  (rad's) 


50  V 


40  V 


30  V 


20  V 


10  V 


140 


40 


120 


160 


Generator  Output  Power  (w) 


21 


Figure  3.5.  Pg-com  characteristic  at  various  values  of  Vtj 
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Figure  3.6.  Pg-Vd  characteristic  at  different  generator  rotational  speed 
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Figure  3.7  shows  the  PMSG  output  power  with  respect  to  duty  ratio  D  at  different 
generator  rotational  speed  when  Vi,us=\50  V.  For  each  rotational  speed,  there  is  only  one 
duty  ratio  that  would  enable  the  system  to  generate  maximum  power.  Therefore,  the 
maximum  power  point  tracking  is  equivalent  to  tracking  the  optimal  duty  ratio.  Control 
of  the  duty  ratio  can  be  achieved  by  a  microcontroller. 
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Figure  3.7.  Pg-D  characteristic  at  different  generator  rotational  speed 


3.4.  PROPOSED  MPPT  CONTROL  ALGORITHM 

Different  hydrokinetic  turbines  have  their  own  Cp-X  characteristic  curve  which 
means  for  different  hydrokinetic  turbines,  the  tip  speed  ratio  for  the  maximum  Cp  is 
different.  It  is  possible  to  measure  the  Cp-X  characteristic  curve  for  each  turbine. 
However,  this  is  not  a  good  solution  for  small-scale  hydrokinetic  turbines  because  the 
cost  to  measure  the  Cp-X  characteristic  curve  for  small-scale  hydrokinetic  turbines  is 
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expensive  compared  to  the  system  cost  itself.  Therefore  a  generic  MPPT  control 
algorithm  is  developed  in  this  thesis  which  automatically  tracks  the  maximum  power 
point  of  the  turbine  without  the  need  to  pre-measure  the  Cp-X  characteristics. 

Figure  3.8  shows  the  flow  chart  of  the  proposed  MPPT  control  algorithm.  P  and  D 
are  the  measured  output  power  and  duty  ratio  of  the  boost  converter,  respectively.  K  is 
the  step  change  in  duty  ratio.  Right  is  a  Boolean  variable  indicating  the  location  of  the 
operating  point.  When  Right  is  true,  the  operating  speed  is  higher  than  the  optimal  rotor 
speed.  When  Right  is  false,  the  operating  speed  is  lower  than  the  optimal  rotor  speed.  At 
the  beginning  of  each  interrupt  cycle,  the  current  cycle  Output  Power  P[n]  is  updated. 
Then  based  on  the  state  of  the  Boolean  variable  Right,  duty  ratio  D  is  changed 
differently.  When  Right  is  true,  if  P[n]>P[n-l],  where  P[n-1]  is  the  output  power  of  the 
previous  cycle,  then  current  cycle  duty  ratio  D[n]=D[n-l]-K  where  D[n-1]  is  the  duty 
ratio  of  the  previous  cycle,  otherwise  if  P[n]<P[n-l],  D[n]=D[n-l]+K  and  Right  needs  to 
be  set  to  false.  When  Right  is  false,  if  P[n]>P[n-l],  D[n]=D[n-l]+K,  otherwise  if 
P[n]<P[n-l],  D[n]=D[n-l]-K  and  Right  needs  to  be  set  to  true.  Since  the  operating  speed 
range  for  hydrokinetic  turbine  is  narrower  than  the  operating  speed  range  for  wind 
turbine  [16],  the  hydrokinetic  turbine  speed  change  cannot  be  too  big  to  cause  a  dramatic 
output  power  change  which  means  the  duty  ratio  step  change,  K,  cannot  be  too  big. 
However,  K  needs  to  be  big  enough  to  cause  a  detectable  output  power  change  for  the 
proposed  MPPT  algorithm  to  work.  Therefore,  the  step  change,  K,  needs  to  be  properly 
tuned  to  each  hydrokinetic  system  for  the  MPPT  control  to  achieve  its  best  performance. 
For  this  project,  K  is  0.025  in  simulation  and  0.0133  in  experimental  test.  In  addition, 
since  the  mechanical  inertia  is  much  higher  than  the  electrical  inertia  of  the  system,  the 
frequency  for  which  P[n]  is  updated  cannot  be  too  high.  This  is  because  P[n]  is  directly 
determined  by  the  steady  state  value  of  com  and  the  mechanical  speed  com  needs  time  to 
reach  its  steady-state  value  because  of  the  high  mechanical  inertia.  Therefore,  the 
frequency  for  which  P[n]  is  updated  needs  to  be  properly  adjusted  as  well  to  each 
hydrokinetic  system  for  the  MPPT  control  to  achieve  its  best  performance.  For  this 
project,  this  frequency  is  10  Hz  simulation  and  experimental  test.  This  MPPT  control 
method  is  similar  to  hill-climbing  used  as  MPPT  control  in  photovoltaic  systems  [23]. 
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Figure  3.8.  Flow  chart  of  the  proposed  MPPT  control  algorithm 


3.5.  SIMULATION  RESULTS 

The  proposed  MPPT  control  is  validated  by  simulation  results  done  in 
Matlab/Simulink  R2010a  and  PLECS.  Figure  3.9  shows  the  block  diagram  of  the 
hydrokinetic  system  built  in  Matlab/Simulink  which  consists  of  the  mechanical  part 
(hydrokinetic  turbine),  the  electrical  part  (generator,  diode  rectifier,  boost  converter  and 
battery)  and  the  control  part  (MPPT). 

Figure  3.10  shows  the  hydrokinetic  turbine  model  built  in  Matlab/Simulink. 
Figure  3.11  shows  the  generator,  diode  rectifier,  boost  converter  and  battery  model  built 
in  PLECS.  The  battery  bank  on  the  load  side  is  simulated  by  a  constant  dc  voltage  source. 
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Figure  3.12  shows  the  block  diagram  of  the  MPPT  control  algorithm  built  in 
Matlab/Simulink. 


Generator  and  Power  Converter 


MPPT  Control 


Figure  3.9.  Block  diagram  of  the  hydrokinetic  system  built  in  Matlab/Simulink 


Figure  3.10.  Block  diagram  of  the  hydrokinetic  turbine  model  built  in  Matlab/Simulink 
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Figure  3.11.  Schematic  of  the  generator  and  power  converters  built  in  PLECS 


Step2 


Figure  3.12.  Block  diagram  of  the  MPPT  control  built  in  Matlab/Simulink 


The  following  four  different  scenarios  are  simulated: 

Case  1  -  Hydrokinetic  system  with  fixed  duty  ratio  and  fixed  water  speed 
Case  2  -  Hydrokinetic  system  with  MPPT  control  and  fixed  water  speed 
Case3  -  Hydrokinetic  system  with  fixed  duty  ratio  and  water  speed  variations 
Case  4  -  Hydrokinetic  system  with  MPPT  control  and  water  speed  variations 
3.5.1.  CASE  1:  Fixed  Water  Speed  and  Duty  Ratio.  The  simplest  scenario  to 
operate  the  hydrokinetic  turbine  system  is  when  there  is  no  input  variations.  Therefore  the 
water  speed  as  well  as  the  duty  ratio  for  the  boost  converter  is  kept  constant.  This 
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operation  mode  does  not  take  advantage  of  the  variable  speed  capability  for  maximum 
power  extraction  of  the  hydrokinetic  turbine.  Figure  3.13  shows  the  simulation  results  of 
water  speed,  turbine  rotational  speed,  Cp,  torque,  duty  ratio  and  the  output  power.  Table 
3.1  shows  the  steady  state  values  of  all  the  parameters  in  Fig.  3.13. 

3.5.2.  CASE  2:  Fixed  Water  Speed  and  MPPT  Control.  The  proposed  MPPT 
control  is  tested  when  there  is  no  water  speed  variations.  This  scenario  is  common  since 
water  speed  tends  to  be  constant.  In  this  case,  the  proposed  MPPT  algorithm 
automatically  tracks  the  optimal  duty  ratio  to  change  the  turbine  rotational  speed  to  its 
optimal  value  where  the  maximum  Cp  occurs.  Figure  3.14  shows  the  simulation  results  of 
water  speed,  turbine  rotational  speed,  Cp,  torque,  duty  ratio  and  the  output  power  under 
this  scenario.  Table  3.2  shows  the  optimal  values  of  all  the  parameters  in  Fig.  3.14.  From 
the  steady-state  values  in  table  3.1  and  the  optimal  values  from  table  3.2,  it  is  proven  that 
at  fixed  water  speed,  hydrokinetic  turbine  system  can  harvest  more  potential  energy  with 
the  proposed  MPPT  control. 

3.5.3.  CASE  3:  Varied  Water  Speed  and  Fixed  Duty  Ratio.  The  third 
operation  scenario  for  the  hydrokinetic  system  is  to  run  it  in  a  varied  water  speed 
environment  with  fixed  duty  ratio.  The  rotational  speed  of  the  turbine  will  change 
because  of  the  water  speed  change.  However,  this  is  not  due  to  the  control  on  the 
electrical  side  of  the  system.  The  water  speed  starts  at  0.9  m/s.  Then  it  is  decreased  by  0.1 
m/s  every  second.  Figure  3.15  shows  the  simulation  results  of  water  speed,  turbine 
rotational  speed,  Cp,  torque,  duty  ratio  and  the  output  power.  Table  3.3  shows  the  steady 
state  values  of  all  the  parameters  in  Fig.  3.15. 

3.5.4.  CASE  4:  Varied  Water  Speed  and  MPPT  Control.  The  most 
complicated  operating  mode  for  the  hydrokinetic  turbine  system  is  to  run  it  in  a  varied 
water  speed  environment  and  in  the  mean  time,  track  the  maximum  power  point  of  the 
hydrokinetic  turbine  under  different  water  speed.  Figure  3.16  shows  the  simulation 
results  of  water  speed,  turbine  rotational  speed,  Cp,  torque,  duty  ratio  and  the  output 
power  under  this  scenario.  Table  3.2  shows  the  optimal  values  of  all  the  parameters  in 
Fig.  3.16.  The  water  speed  starts  at  0.9  m/s  and  the  MPPT  control  changes  the  duty  ratio 
to  its  optimal  value  to  achieve  the  maximum  Cp.  However,  once  the  water  speed 
decreased  to  0.8  m/s,  the  previous  optimal  duty  ratio  is  no  longer  optimal  in  this 
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condition.  The  proposed  MPPT  control  algorithm  can  automatically  adjust  itself  to  the 
change  of  water  speed  and  change  the  value  of  D  until  it  is  optimal  again.  From  the 
steady-state  values  in  table  3.3  and  the  optimal  values  from  table  3.4,  one  conclusion  can 
be  drawn  that  even  with  water  speed  variations  hydrokinetic  turbine  system  can  harvest 
more  potential  energy  with  the  proposed  MPPT  control. 
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Figure  3.13.  Simulation  results  with  fixed  water  speed  and  duty  ratio 
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Figure  3.14.  Simulation  results  with  fixed  water  speed  and  MPPT  control 
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Figure  3.15.  Simulation  results  with  varied  water  speed  and  fixed  duty  ratio 
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Figure  3.16.  Simulation  results  with  varied  water  speed  and  MPPT  control 


33 


Table  3.1.  Steady-state  value  of  the  parameters  of  the  hydrokinetic  system  in  case  1 


Water  speed 

0.9  m/s 

Turbine  rotational  speed 

20.3  rad/s 

Power  coefficient  Cp 

0.38 

Torque 

0.716  N.m 

Duty  ratio 

0.75 

Output  power 

8.96  W 

Table  3.2.  Optimal  value  of  the  parameters  of  the  hydrokinetic  system  in  case  2 


Water  speed 

0.9  m/s 

Optimal  turbine  rotational  speed 

34.45  rad/s 

Maximum  power  coefficient  Cp 

0.395 

Optimal  torque 

0.295  N.m 

Optimal  duty  ratio 

0.55 

Maximum  output  power 

9.84  W 

Table  3.3.  Steady-state  value  of  the  parameters  of  the  hydrokinetic  system  in  case  3 


Water  speed 

0.9  m/s 

0.8  m/s 

0.7  m/s 

0.6  m/s 

Turbine  rotational  speed 

20.3  rad/s 

19.8  rad/s 

19.4  rad/s 

19.1  rad/s 

Power  coefficient  Cp 

0.38 

0.386 

0.391 

0.394 

Torque 

0.716  N.m 

0.352  N.m 

0.243  N.m 

0.157  N.m 

Duty  ratio 

0.75 

0.75 

0.75 

0.75 

Output  power 

8.96  W 

6.53  W 

4.52  W 

2.92  W 
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Table  3.4.  Optimal  value  of  the  parameters  of  the  hydrokinetic  system  in  case  4 


Water  speed 

0.9  m/s 

0.8  m/s 

0.7  m/s 

0.6  m/s 

Optimal  turbine  rotational  speed 

34.45  rad/s 

30.51  rad/s 

24.79  rad/s 

20.9  rad/s 

Maximum  power  coefficient  Cp 

0.395 

0.395 

0.395 

0.395 

Optimal  torque 

0.295  N.m 

0.234  N.m 

0.193  N.m 

0.144  N.m 

Optimal  duty  ratio 

0.55 

0.6 

0.675 

0.725 

Maximum  output  power 

9.84  W 

6.93  W 

4.65  W 

2.94  W 
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3.6.  EXPERIMENTAL  RESULTS 

Experimental  data  are  obtained  through  lab  testing  of  a  hydrokinetic  turbine 
system  shown  in  Fig.  1.3.  The  experimental  setup  is  shown  in  Fig.  3.17.  Experiments  are 
conducted  in  a  controlled  environment  where  the  water  speed  can  be  controlled  and 
measured.  The  maximum  water  speed  for  the  water  tunnel  used  in  this  experimentation  is 
0.9  m/s.  The  constructed  PCB  board  of  the  diode  rectifier  and  the  boost  converter  is 
shown  in  Fig.  3.18.  The  hydrokinetic  turbine  is  used  to  drive  the  PMSG  to  generate 
electrical  power.  However,  the  turbine  rotational  speed  is  too  slow  to  drive  the  PMSG 
used  for  the  experimentation.  Therefore,  a  1:10  speed  increaser  is  used  in  between  the 
hydrokinetic  turbine  and  the  PMSG.  The  output  of  the  PMSG  is  connected  to  the  input  of 
the  diode  rectifier  and  the  output  of  the  boost  converter  is  connected  to  a  constant  voltage 
electronic  load  which  is  used  to  simulate  a  battery  load.  In  order  to  measure  the  input 
voltage  and  current  of  the  boost  converter,  a  Texas  Instruments  TMS320F28335 
microcontroller  is  used.  However,  the  input  range  of  voltage  signals  for  the  ADC 
converter  of  DSP  28335  is  0-3  V.  Therefore  a  circuit  that  changes  the  voltage  and  current 
signal  into  the  allowable  0-3  V  is  required.  This  circuit  is  also  constructed  on  the  PCB 
board  shown  in  Fig.  3.18.  In  addition,  the  rotational  speed  is  measured  by  a  400  ppr 
encoder  that  is  attached  to  the  end  of  the  generator. 

In  order  to  validate  the  proposed  MPPT  control  for  the  hydrokinetic  turbine 
system,  the  following  experiment  procedures  are  adopted.  First  the  P-com  characteristic  of 
the  hydrokinetic  turbine  for  different  water  speed  is  measured  by  manually  varying  the 
duty  ratio  from  0  to  1.  Then  the  MPPT  control  algorithm  is  implemented  in  the  DSP 
28335  microcontroller  to  automatically  track  the  maximum  power  points.  The  maximum 
power  obtained  by  manually  changing  the  duty  ratio  is  then  compared  with  the  maximum 
power  obtained  by  implementing  the  MPPT  control. 

First  the  turbine  power  characteristic  is  tested  when  the  water  speed  is  22  inches/s 
with  the  output  voltage  regulated  at  9.97  V.  The  test  results  are  shown  in  Table  3.5.  Then 
the  turbine  power  characteristic  is  tested  when  the  water  speed  is  26  inches/s  with  the 
same  output  voltage  of  9.97  V  and  the  test  results  are  shown  in  Table  3.6.  Based  on  the 
test  results,  the  turbine  P-co,„  characteristics  under  these  two  different  water  speed  are 
drawn  in  Fig.  3.19.  Also  the  test  results  show  that  when  Ywater=22  inches/s,  the  maximum 
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output  power  is  1.345  W  and  this  happens  when  the  turbine  rotational  speed  is  around  45 
rpm  and  the  optimal  duty  ratio  is  0.467.  When  VWater=26  inches/s,  the  maximum  output 
power  is  2.765  W  and  this  happens  when  the  turbine  rotational  speed  is  around  62  rpm 
and  the  optimal  duty  ratio  is  0.507. 


Table  3.5.  Test  results  of  the  hydrokinetic  system  with  Vwater=22  inches/s,  Vout=9.97V 


Turbine  rotation 

speed  (rpm) 

Duty  ratio 

Output 

Voltage  (V) 

Output 

Current  (A) 

Output 

Power  (W) 

59 

0 

9.97 

0.076 

0.75772 

57 

0.06666667 

9.97 

0.097 

0.96709 

55 

0.13333333 

9.97 

0.108 

1.07676 

53 

0.2 

9.97 

0.119 

1.18643 

51 

0.26666667 

9.97 

0.125 

1.24625 

49 

0.33333333 

9.97 

0.128 

1.27616 

47 

0.4 

9.97 

0.133 

1.32601 

45 

0.46666667 

9.97 

0.135 

1.34595 

43 

0.53333333 

9.97 

0.127 

1.26619 

41 

0.6 

9.97 

0.115 

1.14655 

39 

0.66666667 

9.97 

0.097 

0.96709 

37 

0.73333333 

9.97 

0.083 

0.82751 

34 

0.8 

9.97 

0.061 

0.60817 

32 

0.86666667 

9.97 

0.046 

0.45862 
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Figure  3.17.  Lab  setup  for  experimentation 


Figure  3.18.  PCB  board  of  the  diode  rectifier  and  boost  converter 
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Table  3.6.  Test  results  of  the  hydrokinetic  system  with  VWater=26  inches/s,  Vout=9.97V 


Turbine  rotation 

speed  (rpm) 

Duty  ratio 

Output 

Voltage  (V) 

Output 

Current  (A) 

Output 

Power  (W) 

81 

0 

9.97 

0.0540903 

0.53928 

78 

0.06666667 

9.97 

0.106678 

1.06358 

76 

0.13333333 

9.97 

0.1382307 

1.37816 

74 

0.2 

9.97 

0.1818034 

1.81258 

72 

0.26666667 

9.97 

0.2223711 

2.21704 

69 

0.33333333 

9.97 

0.2464112 

2.45672 

67 

0.4 

9.97 

0.2614363 

2.60652 

64 

0.46666667 

9.97 

0.2764614 

2.75632 

62 

0.507 

9.97 

0.277332 

2.765 

61 

0.53333333 

9.97 

0.2689488 

2.68142 

58 

0.6 

9.97 

0.2494102 

2.48662 

55 

0.66666667 

9.97 

0.2283811 

2.27696 

52 

0.73333333 

9.97 

0.1908185 

1.90246 

48 

0.8 

9.97 

0.1472457 

1.46804 

44 

0.86666667 

9.97 

0.0991655 

0.98868 
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P-Wm  curve  at  Vwater=26 
inches/s 

-  -P-Wm  curve  at  Vwater=22 

inches/s 

- Maximum  Power  at 

Vwater=26inches/s 

. Maximum  Power  at 

Vwater=22inches/s 


Turbine  Rotaional  Speed  (rpm) 

Figure  3.19.  Experimental  results  of  P-com  characteristics  of  hydrokinetic  turbine  under 

two  different  water  speed 


Then  the  hydrokinetic  turbine  system  is  tested  with  the  proposed  MPPT  control 
implemented  in  the  boost  converter  under  Vwater=22  inches/s  and  Vwater=26  inches/s  with 
the  output  voltage  regulated  at  9.97  V.  For  each  water  speed,  10  sample  points  of  the 
output  power  are  collected.  The  test  results  are  shown  in  Table  3.7  and  Table  3.8.  Figure 
3.20  shows  the  test  results  of  the  output  power  with  MPPT  control  compared  with  the 
maximum  output  power  obtained  by  manually  changing  the  duty  ratio  from  0  to  1.  The 
experiment  results  show  that  under  Vwater=22  inches/s,  the  output  power  of  the  system 
oscillates  around  1.345  W  and  that  under  VWater=26  inches/s,  the  output  power  of  the 
system  oscillates  around  2.765  W,  both  of  which  are  the  maximum  output  power  under 
each  water  velocity.  These  experimental  results  validate  the  performance  of  the  proposed 
MPPT  control  for  hydrokinetic  turbine  system  and  its  power  control  capability. 
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Table  3.7.  Test  results  of  the  hydrokinetic  system  with  MPPT  when  Vwater=22  inches/s, 

Vout=9.97V 


sample 

Rotating  speed  (rad/s) 

Output  power  (W) 

1 

45.3 

1.291 

2 

45.9 

1.345 

3 

45.5 

1.325 

4 

46.5 

1.333 

5 

46.3 

1.339 

6 

46.8 

1.314 

7 

45.8 

1.334 

8 

45.5 

1.311 

9 

46.1 

1.344 

10 

46.6 

1.326 

Table  3.8.  Test  results  of  the  hydrokinetic  system  with  MPPT  when  Vwater=26  inches/s, 

Vout=9.97V 


sample 

Rotating  speed  (rad/s) 

Output  power  (W) 

1 

62.5 

2.71 

2 

62.2 

2.754 

3 

62.7 

2.701 

4 

61.9 

2.758 

5 

61.5 

2.712 

6 

61.3 

2.698 

7 

61.8 

2.751 

8 

62 

2.762 

9 

62.4 

2.743 

10 

62.8 

2.733 

Output  Power(W) 
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Sample  points 


♦  Output  Power  with  MPPT 
Control  at  Vwater=26 
inches/s 

■  Output  Power  with  MPPT 
Control  at  Vwater=22 
inches/s 

-  »  Maximum  Output  Power  at 

Vwater=22  inches/s 


Figure  3.20.  Experimental  results  of  output  power  of  hydrokinetic  turbine  with  MPPT 

control  under  two  different  water  speed 
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4.  HIGH-THRUST  LOW-SPEED  PERMANENT  MAGNET 
SYNCHRONOUS  GENERATOR  DESIGN 

4.1.  STRUCTURE  DESIGN  OF  PMSG 

PMSG  has  been  widely  used  in  small-scale  wind  and  hydrokinetic  turbine 
systems  because  of  its  high  reliability,  simple  structure,  low  noise  and  high  power 
density.  In  addition,  due  to  recent  developments  in  permanent  magnet  materials, 
especially  Nd-Fe-B,  high  efficiency  PMSGs  can  be  manufactured  for  turbine  system 
applications  [24].  Traditionally,  the  PMSG  and  the  turbine  are  connected  by  a  gearbox  to 
increase  the  rotational  speed  for  the  PMSG  because  the  rated  operating  speed  of  the 
PMSG  is  usually  much  higher  than  the  operating  speed  of  the  turbines.  However,  it  is 
difficult  to  make  the  efficiency  of  speed  increaser  gearbox  very  high,  which  adds 
additional  power  transmission  losses  in  the  turbine  system.  Therefore,  designing  a  direct 
drive  PMSG  for  hydrokinetic  system  becomes  an  important  issue  in  hydrokinetic  turbine 
systems  [25]. 

In  terms  of  designing  a  direct  drive  PMSG  for  hydrokinetic  turbine  system,  there 
are  a  few  objectives  and  requirements  that  need  to  be  met.  First  of  all,  the  PMSG 
designed  has  to  be  capable  of  producing  the  rated  power  for  the  hydrokinetic  turbine  at 
the  rated  speed  of  the  hydrokinetic  turbine,  which  is  usually  very  low.  This  requires  a 
high-thrust  low-speed  PMSG  design.  Second  of  all,  the  efficiency  of  the  designed  PMSG 
has  to  be  high  enough  so  that  the  direct  drive  configuration  can  surpass  the  gearbox 
configuration  in  terms  of  overall  system  efficiency. 

There  are  two  primary  parts  for  a  PMSG,  the  rotor  which  is  the  rotating  part,  and 
the  stator  which  is  the  stationary  part.  The  main  component  of  the  stator  is  the  stator 
winding  and  the  main  component  of  the  rotor  is  the  permanent  magnet.  Depending  on  the 
placement  of  the  stator  and  the  rotor,  there  are  two  different  structure  designs  for  PMSG. 
Configuration  I  is  the  inner  rotor  and  outer  stator  design  as  shown  in  Fig.  4.1  which  is  a 
4-pole  12-slot  machine.  Configuration  II  is  the  outer  rotor  and  inner  stator  design  as 
shown  in  Fig.  4.2  which  is  a  4-pole  24-slot  machine. 
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Figure  4.1.  Configuration  I:  Inner  rotor  and  outer  stator  design 


Figure  4.2.  Configuration  II:  Outer  rotor  and  inner  stator  design 


Most  of  the  PMSGs  in  the  market  use  configuration  I,  because  when  the  PMSG 
operates,  the  stator  current  going  through  the  stator  winding  would  cause  the  stator 
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temperature  to  go  up.  Therefore,  in  order  to  dissipate  all  the  heat  caused  by  the  stator 
current  more  effectively,  the  stator  is  placed  on  the  outer  side  because  the  rotor  magnets 
act  as  insulators  when  placed  on  the  outer  side.  Furthermore,  in  terms  of  assembling  the 
generator  together  and  coupling  with  the  turbine  shaft,  configuration  I  is  simpler  than 
Configuration  II. 

However,  for  the  hydrokinetic  turbine  applications,  configuration  II  has  several 
advantages.  First  of  all,  the  large  rotor  moment  of  inertia  can  ensure  smooth  operation 
with  small  torque  variations.  In  addition,  configuration  II  can  minimize  the  size  of  stator, 
reducing  the  copper  needed  for  stator  winding  which  ultimately  reduces  the  stator 
resistance  and  stator  inductance,  improving  the  overall  efficiency  of  the  PMSG. 
Therefore,  configuration  II  is  adopted  for  hydrokinetic  turbine  system  applications  taking 
into  account  of  all  the  advantages  it  has  for  this  specific  application. 

4.2.  DESIGN  PROCEDURE  OF  PMSG 

There  are  a  lot  of  different  parameters  that  need  to  be  considered  in  terms  of 
designing  a  PMSG  such  as  the  number  of  poles,  rotor  position,  stator  length  and 
diameters,  rotor  length  and  diameters,  stator  slot  type,  stacking  factor  and  so  on. 
Rotational  Machine  Expert  (RMxprt)  is  a  template-based  electrical  machine  design  tool 
that  provides  fast,  analytical  calculations  of  machine  performance  and  2-D  and  3-D 
geometry  creation  for  detailed  finite  element  calculations  in  ANSYS  Maxwell.  This 
software  streamlines  the  design  procedure  of  PMSG,  therefore  is  used  to  design  a  PMSG 
for  the  hydrokinetic  turbine  system. 

Before  designing  a  PMSG,  some  turbine  information  has  to  be  known  such  as  the 
turbine  maximum  output  power  and  turbine  maximum  rotational  speed.  The  design 
objective  of  the  PMSG  is  to  make  sure  that  the  PMSG  can  generate  the  maximum  output 
power  of  the  turbine  below  the  maximum  rotational  speed  of  the  turbine  with  a  high 
efficiency.  Based  on  the  hydrodynamic  analysis  of  the  hydrokinetic  turbine  that  will  be 
used  in  the  water  tunnel  testing,  the  estimated  maximum  turbine  output  power  is  7.7  w 
and  the  estimated  maximum  turbine  speed  is  240  rpm.  Therefore,  a  PMSG  that  generates 
7.7  w  of  power  below  240  rpm  needs  to  be  designed  which  is  a  high-thrust  low-speed 
machine.  Traditionally,  a  disc  type  machine  (the  diameter  of  the  machine  is  much  larger 
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than  the  length  of  the  machine)  can  get  better  high-thrust  low-speed  performance. 
However,  for  hydrokinetic  turbine  application,  the  size  of  the  PMSG  is  limited  by  the  size 
of  the  hydrokinetic  turbine.  Generally  speaking,  the  diameter  of  the  generator  needs  to  be 
less  than  one  tenth  of  the  diameter  of  the  hub  of  the  hydrokinetic  turbine.  Since  the 
diameter  of  the  hydrokinetic  turbine  that  can  be  tested  for  the  water  tunnel  cannot  exceed 
304.8  mm,  the  maximum  diameter  for  the  PMSG  is  30.48  mm.  However,  this  diameter 
limitation  poses  a  big  challenge  to  design  a  PMSG  to  generate  7.7  w  below  240  rpm.  An 
exception  has  to  be  made  to  increase  the  diameter  of  the  PMSG  which  is  70  mm  in  this 
case.  Figure  4.3  shows  the  designing  template  of  a  PMSG  in  RMxprt.  The  property 
windows  are  accessed  by  clicking  each  of  the  machine  elements,  for  example,  stator, 
rotor  and  shaft.  All  the  parameters  of  the  PMSG  need  to  be  defined  in  the  property 
windows,  and  then  the  designed  machine  needs  to  be  simulated  in  RMxprt.  If  the 
simulation  results  do  not  meet  the  designing  requirements,  then  the  designer  needs  to 
keep  changing  some  of  the  parameters  of  the  machine  and  simulating  the  designed 
machine  until  the  designing  requirements  are  met. 


□  1^3  Brushless  DC  Motor 

nffi  RMxprt  Design  1  (Brushless  Permanent-Magnet  DC  Motor) 
!  i...  Machine 

j--  ©  Circuit 
(■]•••£}  Stator 

MM  i •fi s|ot 

f(J)  Winding 
H  «8>  Rotor 
CO  Pole 
Shaft 

i  0  Analysis 

]•••  0  Setup  1 
I—  0  Setup2 
l 0  Setup3 
Optimetrics 

i  I . Results 

Figure  4.3.  PMSG  design  template  in  RMxprt 
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In  the  Machine  properties  window,  the  machine  type  is  defined  as  Brushless 
Permanent  Magnet  DC  Motor.  The  number  of  poles  is  defined  as  8.  Rotor  position  is 
chosen  to  be  outer  rotor  structure,  the  reason  of  which  is  explained  in  section  4.1.  The 
frictional  loss  and  windage  loss  is  neglected  because  of  their  small  values  and  the 
reference  speed  is  chosen  to  be  the  maximum  rotational  speed  of  the  PMSG  which  is  240 
rpm. 

In  the  rotor  properties  window,  the  outer  diameter  is  defined  as  70  mm.  The  inner 
diameter  is  defined  as  55  mm.  The  length  is  defined  as  76.2  mm  which  is  close  to  the 
diameter  of  the  machine.  The  stacking  factor  is  estimated  to  be  0.95  and  the  pole  type  is 
defined  as  surface  mount  magnet. 

In  the  pole  properties  window,  the  pole  embrace  is  estimated  to  be  0.56,  and  the 
magnet  thickness  is  defined  as  3.175  mm.  The  magnet  material  is  chosen  to  be  Nd-Fe-B. 

In  the  stator  properties  window,  the  outer  diameter  is  defined  as  50  mm  since  it 
has  to  be  smaller  than  the  inner  diameter  of  the  rotor  and  this  gives  an  air  gap  of  2.5  mm. 
The  inner  diameter  is  defined  as  19  mm.  The  length  is  defined  as  76.2  mm.  The  stacking 
factor  is  again  estimated  to  be  0.95.  Total  number  of  slots  is  chosen  to  be  24. 

In  the  slot  properties  window,  the  exact  shape  and  dimension  of  the  slot  can  be 
defined. 

In  the  winding  properties  window,  the  winding  type  is  defined  as  whole-coiled. 
The  number  of  conductors  per  slot  is  chosen  to  be  35. 

Figure  4.4  shows  the  cross  section  of  the  designed  PMSG  after  all  the  parameters 
are  defined.  Figure  4.5  shows  the  stator  winding  configurations  of  the  designed  PMSG. 

4.3.  SIMULATION  RESULTS 

4.3.1.  RMxprt  Simulation.  Once  all  the  parameters  of  the  PMSG  are  defined 
a  solution  setup  needs  to  be  added  to  simulate  the  designed  PMSG.  Solution  setup  is 
defined  in  the  analysis  properties  window.  The  load  type  is  defined  as  constant  speed 
load.  The  rated  output  power  and  speed  is  defined  as  7.7  W  and  240  rpm  respectively. 
The  rated  DC  voltage  is  18  Volts.  Figure  4.6  shows  the  simulation  results  of  the  output 
power  of  the  PMSG  vs.  the  rotational  speed.  The  maximum  output  power  is  1 1 .8  W  at 
145  rpm.  There  are  two  different  operating  points  that  the  generator  can  output  7.7  W, 
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230  rpm  and  60  rpm.  This  simulation  result  meets  the  requirement  of  being  able  to 
generate  7.7  W  below  240  rpm.  Figure  4.7  shows  the  efficiency  vs.  the  rotational  speed. 
At  240  rpm,  the  generator  efficiency  is  around  80%.  Figure  4.8  shows  the  air  gap  flux 
density.  The  peak  value  of  the  air  gap  flux  density  is  650  mT.  Figure  4.9  shows  the 
induced  phase  voltage  and  the  induced  line  voltage  of  the  generator.  The  peak  value  of 
the  phase  voltage  and  the  line  voltage  is  9.3  V  and  16.1  V  respectively.  Figure  4.10 
shows  the  cogging  torque  of  the  designed  PMSG.  Table  4.1  shows  the  simulation  results 
of  some  of  the  important  parameters  of  the  designed  PMSG. 


Table  4.1.  Parameters  of  the  designed  PMSG 


Number  of  poles 

8 

Armature  leakage  inductance 

2.17  mH 

Armature  phase  resistance 

3.32  Ohm 

D-axis  inductance 

3.529  mH 

Q-axis  inductance 

3.529  mH 

Maximum  speed 

300  rpm 

Cogging  torque 

0.0846  N.m 
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\  \  A  1 

IB  1  1 

Figure  4.5.  Winding  layout  of  the  designed  PMSG 
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XY  Plot  1  RMxprtDesign7  A*F 


RSpeed  [rpm] 


Figure  4.6.  Output  power  vs.  rotational  speed 


XY  Plot  2  RMxprtDesign7  AN^F 


RSpeed  [rpm] 

Figure  4.7.  Efficiency  vs.  rotational  speed 


Y1  [V]  i  AirGapFluxDensity  [mTesla] 


50 


XY  Plot  6  RMxprtDesign7  AN^F 


Electrical  Degree  [deg] 
Figure  4.8.  Air  gap  flux  density 
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Figure  4.9.  Induced  phase  and  line  voltage 
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XY  Plot  1  1  RMxprtDesign7  A*F 


The  simulation  results  can  be  verified  by  the  PMSG  design  equations.  First  of  all, 
from  the  pole  pitch  zp  and  the  width  of  the  magnet  tm,  the  magnet  fraction  can  be 
calculated 


a...  =  ■ 


(21) 


Air  gap  flux  density  Bg  can  be  expressed  as 


C,B 

B  = - 


i 


(22) 


In  order  to  calculate  the  air  gap  flux  density  Bg,  flux  concentration  factor  Co, 
permeance  coefficient  Pc,  magnet  leakage  factor  kmi,  carter  coefficient  kc  and  remanence 
B,  needs  to  be  calculated 
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Once  the  air  gap  flux  density  is  obtained,  the  induced  peak  back  EMF  can  be 
expressed  as 

E^=NmBgLRsoNsppnscom  (28) 

where  Nm  is  the  number  of  magnets,  Bg  is  the  air  gap  flux  density,  L  is  the  length 
of  the  stator,  Rso  is  the  stator  outer  radius,  Nspp  is  the  number  of  slots  per  pole  per  phase, 
ns  is  the  number  of  turns  per  slot  and  com  is  the  rotational  speed  of  the  generator. 

After  applying  the  parameters  of  the  designed  PMSG  into  the  equations  above, 
the  peak  back  EMF  is  7.09  V  which  is  close  to  the  simulation  result  9.3  V.  The  difference 
of  these  two  results  may  be  caused  by  the  approximation  used  in  the  equations  above. 
Also,  the  distribution  factor,  pitch  factor  and  the  skew  factor  are  not  taken  into  account 
during  the  above  calculations. 

4.3.2.  Magnetostatic  Simulation  in  Maxwell  2D.  Though  RMxprt  is  able  to 
simulate  the  performance  of  the  designed  PMSG  in  terms  of  output  power,  efficiency, 
induced  voltage,  and  other  dynamic  performance,  it  does  not  give  the  magnetostatic 
properties  of  the  designed  PMSG.  In  order  to  simulate  the  magnetostatic  properties  such 
as  flux  distribution  and  flux  density  distribution  of  the  designed  PMSG  in  RXxprt,  a 
finite  element  analysis  (FEA)  needs  to  be  conducted.  Therefore,  a  Maxwell  2D  model 
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needs  to  be  designed.  Figure  4.11  shows  the  cross  section  view  of  the  designed  PMSG  in 
Maxwell  2D.  The  parameters  and  dimensions  of  this  PMSG  is  the  same  as  the  one 
designed  using  RMxprt.  Figure  4.12  shows  the  flux  distribution  of  the  designed  PMSG  in 
Maxwell  2D.  Figure  4.13  shows  the  flux  density  distribution  of  the  designed  PMSG  in 
Maxwell  2D.  The  maximum  flux  density  is  1.28  T  and  this  high  flux  density  only  exists 
at  the  edge  of  the  magnets  where  the  flux  paths  through  the  outer  steel  part  of  the  rotor. 
Therefore,  this  PMSG  design  does  not  have  flux  saturation  issues. 


Stator  Lamination  (Electrical  Steel) 
Rotor  (Steel) 

Permanent  Magnet  (NdFcB) 


Winding 


(Stainless  Steel) 


Figure  4.11.  Cross  section  view  of  the  designed  PMSG  in  Maxwell  2D 
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Figure  4.12.  Flux  line  distribution  of  the  designed  PMSG  in  Maxwell  2D 
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Figure  4.13.  Flux  density  distribution  of  the  designed  PMSG  in  Maxwell  2D 
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4.4.  TEST  RESULTS 

In  order  to  test  the  performance  of  the  designed  PMSG,  the  physical  motor  is  built 
as  shown  in  Fig.  4.14.  One  of  the  most  challenging  difficulties  encountered  during 
construction  of  the  motor  is  that  the  stator  lamination  pieces  are  very  sharp  and  they  tend 
to  cut  into  the  insulation  layers  of  the  magnetic  wire  causing  multiple  shorts  among  the 
different  phases.  In  order  to  avoid  this  problem,  additional  insulation  needs  to  be 
provided  on  the  windings.  In  this  case,  electrical  tape  is  used  in  the  stator  slots  to  provide 
the  additional  protection  between  the  laminations  and  the  windings.  However  the 
available  space  in  the  slots  for  windings  is  reduced.  When  the  motor  is  finally 
constructed,  there  are  only  10  turns  per  slot  instead  of  the  35  turns  per  slot  as  designed 
previously.  The  actual  phase  inductance  and  phase  resistance  is  measured  which  is  249 
pH  and  1.253  Ohms  respectively. 


Figure  4.14.  PMSG  built  based  on  the  RMxprt  design 


An  induction  motor  is  used  to  drive  the  PMSG  to  test  its  performance.  Figure  4.15 
(a)-Figure  4.22  (a)  shows  the  experiment  results  of  generated  three  phase  voltages  under 
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different  rotational  speed.  Figure  4.15  (b)-Figure  4.22  (b)  shows  the  simulation  results  of 
induced  phase  voltages  under  the  corresponding  rotational  speed  in  experiment. 


Figure  4.15  (a)  Experimental  result  of  three  phase  voltage  at  84  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.15  (b)  Simulation  result  of  induced  phase  voltage  at  84  rpm 


InducedPhaseVoltage  [V] 
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Figure  4.16  (a)  Experimental  result  of  three  phase  voltage  at  171  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.16  (b)  Simulation  result  of  three  phase  voltage  at  171  rpm 


InducedPhaseVoltage  [V] 
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Figure  4.17  (a)  Experimental  result  of  three  phase  voltage  at  258  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.17  (b)  Simulation  result  of  three  phase  voltage  at  258  rpm 


InducedPhaseVoltage  [V] 
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Figure  4.18  (a)  Experimental  result  of  three  phase  voltage  at  345  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.18  (b)  Simulation  result  of  three  phase  voltage  at  345  rpm 


InducedPhaseVoltage  [V] 


60 


2B  2  00  V 

m 

« 

m  2.00  V 

i 

IM  lOOmV 

4  QOliS 

9  1  rfqueivy 
9  CMS 

2-  i:  f  •  . 

jfcnvis 

'  rrqur  ncv 

.fas 

V»|UP 

;k  r.s  II? 

)  154  V 

Moon 

n  4% 

4  206 

Min 

29  69 

.1  IS4 

Mix 
■*0  1  1 

5  012 

Vfd  Wv 

7  224 

771  Sni 

2K.4B  H/ 
3.032  V 

39.44 

4  102 

21  4B 

3.032 

4.92S 

7.210 

KOI  Ont 

j  10  Oim 

ounoo 

10  .OMS/s  1  1 

1  IM  piunts  1 

Cuupliwj 

m 


Tcrmifutlon 

by 

AC  TPF0S00 


inwfi 

On  (on) 


BincJwKUh 

lull 


2  Ubtl 


J  MOV 


6  VUf  201 2 
01  19 


Figure  4.19  (a)  Experimental  result  of  three  phase  voltage  at  432  rpm 
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Figure  4.19  (b)  Simulation  result  of  three  phase  voltage  at  432  rpm 


InducedPhaseVoltage  [V] 


61 


•  2  0#  V 

*  - 

A  2.00  V 

igji  lOOmV 

4  OOMl 

Vali* 

Moan 

Min 

Max 

Std  I>rv  i 

O  1 cwiijb  ncv 

34.23  M/ 

34  23 

34  23 

34  23 

o.ooo  L 

•  i:m> 

3.733  V 

3.723 

3.733 

O.OOG 

I 

fl§  t  y 

34  10  H? 

14  10 

J4  10 

14  It 

0  OHO 

1,07)  V 

1  07) 

1  0?) 

i  071 

0,0110 

0  Mm  joij 

01  no  47 


Figure  4.20  (a)  Experimental  result  of  three  phase  voltage  at  519  rpm 
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Figure  4.20  (b)  Simulation  result  of  three  phase  voltage  at  519  rpm 
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Figure  4.21  (a)  Experimental  result  of  three  phase  voltage  at  606  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.21  (b)  Simulation  result  of  three  phase  voltage  at  606  rpm 


InducedPhaseVoltage  [V] 
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Figure  4.22  (a)  Experimental  result  of  three  phase  voltage  at  693  rpm 


XY  Plot  8  RMxprtDesign7 


Figure  4.22  (b)  Simulation  result  of  three  phase  voltage  at  693  rpm 
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Figure  4.15-Figure  4.22  shows  that  the  experiment  results  of  the  generated  phase 
voltage  is  the  same  as  the  simulation  results  in  RMxprt.  Table  4.2  shows  the  experimental 
results  of  the  rms  value  of  the  induced  phase  voltage,  back  EMF  constant,  electrical 
frequency  coe,  mechanical  frequency  com,  and  the  ratio  of  coe  over  com  under  different 
rotational  speed.  Since  it’s  an  8  pole  machine,  the  ratio  of  co/com  should  be  the  number  of 
pole  pairs  which  is  4.  This  is  validated  by  the  results  shown  in  Table  4.2  where  coe/com  is 
close  to  4  under  different  rotational  speed.  Figure  4.23  shows  the  relationship  between 
the  rms  value  of  the  induced  phase  voltage  and  the  rotational  speed.  Based  on  Figure 
4.23,  the  rms  value  of  the  induced  phase  voltage  and  the  rotational  speed  of  the  PMSG  is 
almost  a  linear  relationship  and  this  relationship  can  be  characterized  by  the  parameter 
called  back  EMF  constant.  According  to  table  4.2,  the  back  EMF  constant  of  this 
machine  is  around  0.007  V/rpm. 


Table  4.2.  Test  results  of  PMSG  under  different  rotational  speed 


Speed  (rpm) 

Phase  voltage  rms  (V) 

BACK  EMF  V/rpm 

COe 

(De/ COm 

84 

0.584415827 

0.006957331 

5.55 

1.4 

3.96905 

171 

1.204479561 

0.00704374 

11.24 

2.85 

3.94503 

258 

1.82442574 

0.007071418 

17.04 

4.3 

3.96357 

345 

2.45248916 

0.007108664 

22.9 

5.75 

3.98319 

432 

3.076042872 

0.00712047 

28.65 

7.2 

3.98009 

519 

3.681418494 

0.007093292 

34.12 

8.65 

3.94528 

606 

4.304191447 

0.007102626 

39.48 

10.1 

3.90957 

693 

4.933487348 

0.007119029 

45.85 

11.55 

3.96999 
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5.  CONCLUSIONS  AND  FUTURE  WORK 


5.1.  CONCLUSION 

River-based  hydrokinetic  turbine  system  is  studied  and  discussed  in  terms  of 
system  level  design  and  power  flow  control  method  in  this  thesis.  Hydrokinetic  turbine 
power  characteristics  is  also  analyzed  and  its  similar  capability  of  variable  speed 
operation  for  maximum  power  extraction  has  led  to  the  proposal  of  maximum  power 
point  tracking  (MPPT)  control  for  hydrokinetic  turbines. 

A  complete  mathematical  model  of  the  hydrokinetic  system  is  developed.  This 
model  analyzes  the  interaction  between  the  mechanical  part  of  the  system  and  the 
electrical  part  of  the  system.  Based  on  the  developed  model  of  the  hydrokinetic  turbine 
system,  a  boost  converter  is  used  to  implement  the  proposed  MPPT  control.  The 
developed  MPPT  control  is  similar  to  the  hill-climbing  method  used  as  MPPT  control  in 
photovoltaic  systems.  In  addition,  the  proposed  MPPT  control  strategy  is  generic  to  all 
the  hydrokinetic  turbine  systems.  Simulation  and  experiment  results  validate  the 
performance  of  the  proposed  control  strategy. 

In  order  to  improve  the  efficiency  of  the  power  transmission  of  the  hydrokinetic 
system,  a  low-speed  high-thrust  permanent  magnet  synchronous  generator  (PMSG)  is 
recommended  for  the  small  scale  stand-alone  hydrokinetic  turbine  systems.  The  adoption 
of  a  low-speed  high-thrust  generator  can  eliminate  the  use  of  a  speed  increaser,  thus 
reducing  the  transmission  power  losses  imposed  by  the  speed  increaser.  A  low-speed 
high-thrust  PMSG  is  also  designed  for  the  hydrokinetic  turbine  system  using  RMxprt.  A 
Maxwell  2D  model  of  the  designed  PMSG  is  also  developed  to  simulate  the 
magnetostatic  performance  of  the  machine.  A  physical  prototype  of  the  designed  PMSG 
is  built  and  tested  as  well. 

5.2.  FUTURE  WORK 

River-based  hydrokinetic  turbine  system  is  still  at  its  elementary  stage  in  terms  of 
commercialization.  There  are  still  a  lot  to  be  done.  In  terms  of  structure,  a  new  design  is 
necessary  for  the  integration  of  the  hydrokinetic  turbine  and  the  electrical  system  such  as 
the  PMSG  and  the  power  converters.  In  terms  of  electrical  part  for  the  system,  grid 
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interface  can  be  developed  so  that  the  extra  amount  of  potential  power  from  the 
hydrokinetic  turbine  can  be  injected  into  the  grid,  saving  the  customers  money  and 
making  the  hydrokinetic  turbine  system  more  cost-effective.  In  addition,  the  system  can 
be  scaled  up  so  that  it  has  more  power  capability. 


APPENDIX  A. 

PRINTED  CIRCUIT  BOARD  DESIGN 
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This  appendix  includes  the  schematic  and  board  layout  of  the  printed  circuit  board 
described  in  Section  3.6.  The  schematic  diagram  is  shown  in  Figure  A.l.  The  physical 
board  layout  pictured  in  Figure  3.17,  is  broken  into  two  parts  to  show  the  top  and  the 
bottom  layers  of  the  board  separately.  The  board  layout  for  the  top  layer  is  shown  in 
Figure  A. 2  and  the  board  layout  for  the  bottom  layer  is  shown  in  Figure  A. 3. 


Figure  A.l.  Printed  circuit  board  schematic 


Figure  A. 2.  Physical  board  layout  of  top  layer 


Figure  A. 3.  Physical  board  layout  of  bottom  layer 


APPENDIX  B. 

DSP  CODE  FOR  THE  PROPOSED  MPPT  CONTROL  ALGORITHM 
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This  appendix  includes  the  C  code  of  the  proposed  MPPT  algorithm  programmed 
into  the  DSP28335  microcontroller.  There  are  detailed  comments  of  explanation  for  each 
block  of  codes.  The  comments  are  denoted  by  text  which  is  preceded  with  two  forward 
slash  symbols  (//). 


//  Copyright  by  Hailong  Zhou 
//  All  Rights  Reserved 

//  MPPT  Control  Algorithm  for  Hydrokinetic  Turbine 
//############################################################################# 
//  Code  Description: 

//  This  code  is  used  to  implement  the  MPPT  control  in  the  boost  converter 

//  for  the  hydrokinetic  turbine  system.  V_decimal  and  I_decimal  is  the 
//  measured  voltage  and  current.  These  two  are  updated  every  0.01  ms.  However, 
//  the  power  Pnow  is  updated  every  0.1  s.  The  duty  ratio  is  changed  by  varing 
//  the  value  of  EPwmlRegs . CMPA.half . CMPA.  The  MPPT  control  and  the  speed 
//  sensing  is  realized  in  the  interrupt  service  of  timerO. 
//############################################################################# 

//############################################################################# 

#include  "DSP2 833x_Device . h" 

//  External  function  prototypes 

extern  void  InitAdc (void) ; 
extern  void  InitSysCtrl (void) ; 
extern  void  InitPieCtrl (void) ; 
extern  void  InitPieVectTable (void) ; 
extern  void  InitCpuTimer s (void) ; 

extern  void  Conf igCpuTimer (struct  CPUTIMER_VARS  *,  float,  float); 

//  Prototype  statements  for  functions  found  within  this  file. 

void  Gpio_select (void) ; 

void  Setup_ePWMlA (void)  ; 

void  Setup_ePWM2 (void) ; 

void  Setup_ADC (void) ; 

interrupt  void  cpu_timer 0_isr (void) ;  //CPU  TIMERO 

interrupt  void  adc^isr (void) ;  //  ADC  End  of  Sequence  ISR 

//Global  variable. 

int  count  =  0; 

int  start-up  =  1; 

int  up  =  0 ; 

float  P  now  =  0 ; 

float  Pbefore  =  0; 

unsigned  int  Voltage  V; 

unsigned  int  Voltagel; 

unsigned  int  Vg; 

float  V_decimal  =  0; 

float  I_decimal  =  0; 

float  P_decimal_now  =  0; 

float  P_decimal_before  =  0; 

float  Duty_ratio  =  0; 

float  Dnumber  =  0; 

float  pi  =  3.141592653; 

float  T1  =  0; 

float  state_now  =  0; 

float  state  before  =  0; 
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float  fe  =  0; 
float  We  =  0; 
float  Werpm  =  0; 

float  Wm  =  0 ; 
float  Wmrpm  =  0; 
float  pulse_count  =  0; 
float  Torque  =  0; 
float  Rps  =  0; 


//############################################################################# 

/ /  main  code 

//############################################################################# 

void  main (void) 

{ 

InitSysCtr 1  ( ) ;  //  Basic  Core  Init  from  DSP2833x_SysCtrl . c 

EALLOW; 


SysCtrlRegs . WDCR= 
EDIS; 

DINT; 

Gpio_select ( )  ; 
Setup_ePWMlA ( )  ; 
InitPieCtr 1 ( )  ; 
InitPieVectTable ( ) 
Setup_ADC()  ; 
Setup_ePWM2 ()  ; 


OxOOAF;  //  Re-enable  the  watchdog 
//  OxOOAF  to  NOT  disable  the  Watchdog,  Prescaler  =  64 

//  Disable  all  interrupts 

//  GPIO  set  up 

//  Init  of  ePWMIA 

//  Basic  setup  of  PIE  table;  from  DSP2833x_PieCtrl . c 
;  //  Default  ISR's  in  PIE 

//  Init  of  ADC 
//  Init  of  ePWM2 


EALLOW; 

PieVectTable . TINT0  =  &cpu_timer0_isr ; 

PieVectTable . ADCINT  =  &adc_isr; 

EDIS; 

InitCpuTimer s ( ) ;  //  Basic  setup  CPU  TimerO,  1  and  2 

Conf igCpuTimer (SCpuTimerO, 150, 10) ;  //  Timer  configuration,  0.01ms 

PieCtrlRegs. PIEIER1 . bit . INTx6  =1;  //  ADC 

PieCtr IRegs . PIEIER1 . bit . INTx7  =1;  //  CPU  Timer 


IER  |=1; 


EINT; 

ERTM  ; 

CpuTimerORegs . TCR. bit . TSS  =  0;  //  Start  timerO 

while ( 1 ) 


EALLOW; 

SysCtrlRegs . WDKEY  =  0x55;  //  Service  WD  #1 

EDIS; 

} 

} 


void  Gpio_select (void)  //  GPIO  set  up 
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EALLOW; 

GpioCtrlRegs . GPAMUX1 . all  =  0; 
Puropse  I/O 

GpioCtrlRegs. GPAMUX1. bit. GPIOO  =  1; 

GpioCtrlRegs .GPAMUX2 .all  =  0; 
Purpose  I/O 

GpioCtrlRegs .GPBMUX1 .all  =  0; 

Purpose  I/O 

GpioCtrlRegs .GPBMUX2 .all  =  0; 
Purpose  I/O 

GpioCtrlRegs .GPCMUX1 .all  =  0; 
Purpose  I/O 

GpioCtrlRegs .GPCMUX2 .all  =  0; 
Purpose  I/O 

GpioCtrlRegs . GPADIR. all  =  0; 
GpioCtrlRegs . GPADIR. bit . GPIOO  = 
GpioCtrlRegs . GPBDIR . all  =  0; 
GpioCtrlRegs . GPCDIR. all  =  0; 

EDIS; 

} 

void  Setup_ePWMlA (void) 

{ 

EPwmlRegs . TBCTL . bit . CLKDIV  =  0; 
EPwmlRegs . TBCTL . bit . HSPCLKDIV  =  0; 
EPwmlRegs. TBCTL. bit. CTRMODE  =  2; 

EPwmlRegs. AQCTLA. all  =  0x0060; 


EPwmlRegs. TBPRD  =  750; 
EPwmlRegs. CMPA.half.CMPA  =  7  50; 

} 

void  Setup_ePWM2 (void) 

{ 

EPwm2Regs . TBCTL. all 
/* 
bit 
bit 
bit 
bit 
bit 
bit 
bit 
bit 
bit 
*/ 


// 

GPI015 

. . .  GPIOO  = 

=  General 

// 

ePWMIA 

active 

// 

GPI031 

. . .  GPI016 

=  General 

// 

GPI047 

. . .  GPI032 

=  General 

// 

GPI063 

. . .  GPI048 

=  General 

// 

GPI079 

. . .  GPI064 

=  General 

// 

GPI087 

. . .  GPIO80 

=  General 

1;  //  GPIOO  as  output  PWM 

//  GPI063-32  as  inputs 
//  GPI087-64  as  inputs 

//Set  up  ePWMIA  100K 

//  CLKDIV  =  1 
//  HSPCLKDIV  =  2 
//up  -  down  mode 

//  Set  ePWMIA  on  CMPA  up 
//  Clear  ePWMIA  on  CMPA  down 

//  lOOKHz  -  PWM  signal 
//  0%  duty  cycle  first 

//Set  up  ePWM2  100K 

//  Configure  timer  control  register 

11  =  ignore  emulation  suspend 
count  down  after  sync  event 


0xC030 ; 


15-14 

11 : 

FREE/ SOFT 

13 

0: 

PHSDIR,  0 

12-10 

000: 

CLKDIV,  0 

9-7 

000: 

HSPCLKDIV 

6 

0: 

SWFSYNC , 

5-4 

11  : 

SYNCOSEL, 

3 

0: 

PRDLD,  0  = 

2 

0: 

PHSEN,  0  ^ 

o 

1 

\ — 1 

00: 

CTRMODE, 

00  =>  TBCLK  =  HSPCLK/1 
,  000  =>  HSPCLK  =  SYSCLKOUT/1 
0  =  no  software  sync  produced 
11  =  sync-out  disabled 
=  reload  PRD  on  counter=0 


00 


count  up  mode 


EPwm2Regs . TBPRD  =  2999;  //  TBPRD  +1  =  TPWM  /  (HSPCLKDIV  *  CLKDIV  * 

TSYSCLK) 

//  =  20  M  /  6.667  ns 

//100K 


EPwm2Regs.ETPS.all  =  0x0100; 

/* 


//  Configure  ADC  start  by  ePWM2 


bit  15-14 

00: 

EPWMxSOCB,  read-only 

bit  13-12 

00: 

SOCBPRD,  don't  care 

bit  11-10 

00: 

EPWMxSOCA,  read-only 

bit  9-8 

01 : 

SOCAPRD,  01  =  generate 

bit  7-4 

0000: 

reserved 
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} 


bit  3-2  00:  INTCNT,  don't  care 

bit  1-0  00:  INTPRD,  don't  care 

*/ 

EPwm2Regs . CMPA. half . CMPA  =  15000;  //  0%  duty  cycle  first 

EPwm2Regs.ETSEL.all  =  OxOAOO;  //  Enable  SOCA  to  ADC 

/* 


bit 

15 

0: 

SOCBEN,  0  =  disable 

SOCB 

bit 

14-12 

000: 

SOCBSEL,  don't  care 

bit 

11 

1 : 

SOCAEN,  1  =  enable 

SOCA 

bit 

10-8 

010: 

SOCASEL,  010  =  SOCA 

on  PRD  event 

bit 

7-4 

0000: 

reserved 

bit 

3 

0: 

INTEN,  0  =  disable 

interrupt 

bit 

o 

1 

C\] 

000: 

INTSEL,  don't  care 

*/ 


void  Setup_ADC (void) 

{ 

InitAdcO;  //  Basic  ADC  setup,  incl .  calibration 


AdcRegs . ADCTRL1 . all 
AdcRegs . ADCTRL1 . bit 
AdcRegs . ADCTRL1 . bit 
AdcRegs . ADCTRL1 . bit 
AdcRegs . ADCTRL1 . bit 

AdcRegs . ADCTRL2 . all 
AdcRegs . ADCTRL2 . bit 
AdcRegs . ADCTRL2 . bit 

trigger 

AdcRegs . ADCTRL2 . bit 
of  sequence 


=  0; 


ACQ  PS  =  7; 

// 

SEQ  CASC  =  1; 

// 

CPS  =  0; 

// 

CONT  RUN  =  0; 

// 

=  0; 

INT  ENA  SEQ1  = 

1; 

// 

EPWM  SOCA  SEQ1 

=  l; 

// 

INT  MOD  SEQ1  = 

0; 

// 

7  =  8  x  ADCCLK 
l=cascaded  sequencer 
divide  by  1 
single  run  mode 


l=enable  SEQ1  interrupt 
1=SEQ1  start  from  ePWMJSOCA 

0=  interrupt  after  every  end 


AdcRegs . ADCTRL3 . bi t . ADCCLKPS 
ADCCLKPS 

DSP2 833x_SysCtrl .  c) 


=  3;  //  ADC  clock:  FCLK  =  HSPCLK  /  2  * 
//  HSPCLK  =  75MHz  (see 
//  FCLK  =12.5  MHz 


AdcRegs. ADCMAXCONV. all  =  3; 

AdcRegs. ADCCHSELSEQ1 .bit. CONV00 
AdcRegs . ADCCHSELSEQ1 . bit . CONV01 
AdcRegs . ADCCHSELSEQ1 . bit . CONV02 

} 

interrupt  void  cpu_timer 0_isr (void) 


//  3  conversions  from  Sequencer  1 


0; 

// 

Setup 

ADCINA0 

as 

1st 

SEQ1 

conv . 

1; 

// 

Setup 

ADCINA1 

as 

2nd 

SEQ1 

conv 

2; 

// 

Setup 

ADCINA2 

as 

3rd 

SEQ1 

conv 

CpuTimerO . InterruptCount++ ; 
count++; 

EALLOW; 

SysCtrlRegs . WDKEY  =  OxAA;  //  Service  WD  #2 
EDIS; 

//########################################################################### 

//  MPPT  control 

//########################################################################### 

if (start-up) 


if (EPwmlRegs .CMPA. half .CMPA  >300) 

{ 


if (CpuTimerO . InterruptCount%500==0) 
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{ 

EPwmlRegs . CMPA . hal f . CMPA  =  EPwmlRegs . CMPA. half . CMPA-1 ; 

} 

} 

else 

{ 

start_up  =  0; 
up  =  1  ; 

} 

} 

else 

{ 

if  (count%10000==0) 

{ 

switch (up) 

{ 

case  1 : 

P_now  =  (Voltage_V/10) * (Voltage_I/10) ; 

P_decimal_now  =  V_decimal* I_decimal ; 
if (P_decimal_now>P_decimal_bef ore) 

{ 

P_decimal_before  =  P_decimal_now; 

EPwmlRegs. CMPA. half .CMPA  =  EPwmlRegs . CMPA. half . CMPA+10 ; 

} 

else 


P_decimal_before  =  P_decimal_now; 

EPwmlRegs. CMPA. half .CMPA  =  EPwmlRegs . CMPA. half . CMPA-10 ; 
up  =  0; 

} 

break; 
case  0 : 

P_decimal_now  =  V_decimal* I_decimal  ; 
if (P_decimal_now>P_decimal_bef ore) 

{ 

P_decimal_before  =  P_decimal_now; 

EPwmlRegs. CMPA. half .CMPA  =  EPwmlRegs . CMPA. half . CMPA-10 ; 

} 

else 


P_decimal_before  =  P_decimal_now; 

EPwmlRegs. CMPA. half .CMPA  =  EPwmlRegs . CMPA. half . CMPA+10 ; 
up  =  1  ; 

} 

break; 

} 

count  =  0; 

} 

} 

//##################################################################### 

//  Speed  sensing  of  the  generator 

//##################################################################### 

state_now  =  GpioDataRegs . GPADAT . bit . GPI09; 


T1  =  Tl+0. 00001; 


if ( ( state_before==0) && (state_now==l) ) 

{ 

pulse_count  =  pulse_count+l ; 
if (pulse_count==1000 ) 

{ 


fe  =  pulse_count/Tl ; 
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Rps  =  fe/400; 

Werpm  =  Rps* 60; 

Wmrpm  =  Werpm/2; 

Wm  =  Wmrpm*2*pi/60; 

Torque  =  P_decimal_now/Wm; 

T1  =  0; 

pulse_count  =  0; 

} 

} 

state_before  =  state_now; 

Dnumber  =  EPwmlRegs . CMPA.half . CMPA; 
Duty_ratio  =  Dnumber/750; 

PieCtr IRegs . PIEACK .all  =  PIEACK_GR0UP1 ; 


interrupt  void  adc_isr (void) 

{ 

Voltage_V  =  AdcMirror .ADCRESULTO; //  Voltage  measurement 

Voltage_I  =  AdcMirror.ADCRESULTl;  //  Current  measurement 

V_decimal  =  Voltage_V*0. 0212688;  //  Voltage  value  in  decimal 

I_decimal  =  Voltage_I*0 . 00012152 ;  //  Current  value  in  decimal 

//  Reinitialize  for  next  ADC  sequence 

AdcRegs . ADCTRL2 . bit . RST_SEQ1  =1;  //  Reset  SEQ1 

AdcRegs.ADCST.bit. INTJ3EQ1JCLR  =1;  //  Clear  INT  SEQ1  bit 

PieCtrlRegs . PIEACK . all  =  PIEACK_GR0UP1 ;  //  Acknowledge  interrupt  to  PIE 

} 

//===================================================================== 

//  End  of  Code. 

//===================================================================== 
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